Condensation and Clustering in the Driven Pair Exclusion Process 
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We investigate particle condensation in a driven pair exclusion process on one- and two- dimen- 
sional lattices under the periodic boundary condition. The model describes a biased hopping of 
particles subject to a pair exclusion constraint that each particle cannot stay at a same site with 
its pre-assigned partner. The pair exclusion causes a mesoscopic condensation characterized by the 
scaling of the condensate size m con ~ N" and the number of condensates N con ~ N a with the total 
number of sites N. Those condensates are distributed randomly without hopping bias. We find 
that the hopping bias generates a spatial correlation among condensates so that a cluster of con- 
densates appears. Especially, the cluster has an anisotropic shape in the two-dimensional system. 
The mesoscopic condensation and the clustering are studied by means of numerical simulations. 
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I. INTRODUCTION 

In a driven diffusive system, hopping bias in dynam- 
ics drives a system out of equilibrium and can lead to a 
variety of effects. The hopping bias is irrelevant for non- 
interacting or single-particle systems. Biased diffusion of 
a particle, for example, is equivalent to unbiased diffusion 
in a co-moving frame. Particle interaction can make a dif- 
ference. Consider, as an example, the asymmetric simple 
exclusion process (ASEP) [l| in one dimension (ID), in 
which particles hop under the constraint that each site 
can be occupied by at most a single particle [exclusion 
interaction]. When particles hop symmetrically, the sys- 
tem belongs to the Edwards- Wilkinson universality class 
in which the dynamic correlation length £ grows in time 
t as £ ~ t 2 Q. With a hopping bias, however, the system 
exhibits a scaling £ ~ i 3 / 2 of the Kardar-Parisi-Zhang 
universality class [3j. 

In general, the hopping bias can also change the nature 
of a stationary state probability distribution. We focus 
on systems displaying a condensation transition. Con- 
densation occurs in various forms such as Bosc-Einstein 
condensation, traffic jams, hub formation in evolving net- 
works, and so on [4|. Condensation can be studied in the 
context of driven diffusive systems [J, Q , and the effect 
of hopping bias has been studied in those systems [6H9| . 

The zero range process (ZRP) is useful in studying 
particle condensation [4|, [lfl, [ll|. Consider, for simplic- 
ity, a one-dimensional ring of N sites i — 1,2, ■•■ ,7V 
under a periodic boundary condition. There are M par- 
ticles, and the occupation number at site i is denoted as 
rrii = 0, 1, 2, • • •. A particle at site i selected randomly 
decreases by one (m, -> m; - 1) at a rate u(rrii) and 
then hops to a neighboring site j = i + 1 with probability 
p or j = i — 1 with probability q = 1 — p. The jump- 



ing rate depends on the occupation number at a source 
site, whose functional form reflects the nature of on-site 
interactions among particles. For example, if particles 
hop independently, the jumping rate should be linearly 
proportional to the occupation number, u(m) ex m. A 
jumping rate function u(m) growing sublinearly or de- 
creasing in m corresponds to an attractive interaction. 

The stationary state of the ZRP is known exactly. The 
probability P st (m) to find the system in a configuration 
m = (mi, ■ ■ ■ , itil) in the stationary state is given by a 
product form 
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with a normalization factor Z(N, M). The exact solution 
allows one to understand the condition for condensation. 
Suppose that the jumping rate is given by u(m) = 1 + 
b/m. When b < 2, particles are distributed uniformly at 
any value of particle density p = M/N. On the other 
hand, when b > 2, there emerges a single macroscopic 
condensate of size m con = (p — p c )N for p > p c with 
p c = 1/(6 — 2). The ZRP has the same stationary state 
at all values of p and q irrespective of the hopping bias. 

An interesting generalization of the ZRP was consid- 
ered in Ref. [6|, where the particle jumping rate is con- 
stant, but a jump is accepted with probability v(m) de- 
pending on the occupation number m at target sites. 
Such a model is called the target process (TP). With- 
out a hopping bias (p = q), the stationary state of the 
TP is mapped to that of the ZRP with jumping rate 
u(m + 1) = l/v(m). Hence, the unbiased TP displays 
macroscopic condensation under an appropriate condi- 
tion, e.g., v(m) = 1/(1 + b/(m + 1)) with 6 > 2 and 
p > 1/(6—2). However, when a hopping bias with 
p =/= q exists, the mapping breaks down, and the sta- 
tionary state of the TP is not given by the product form. 
Moreover, the hopping bias destroys the condensation [fj. 

The conserved mass aggregation (CMA) model also 



exhibits a condensation phenomenon Q. In this model, 
all particles at a site i hop (rrii -> 0) to a neighbor- 
ing site at the unit rate or a single particle is chipped 
off (rrii — > rrii — 1) to a neighboring site at a constant 
rate u. Particles at different sites aggregate through the 
former process while they are scattered out by the lat- 
ter. Competition between them results in condensation 
when p > Pc{^) = \/uT+T — 1 in all dimensions when 
the hopping and the chipping are symmetric 0- How- 
ever, a bias in the hopping and the chipping was shown 
to inhibit condensation in one-dimension [8[ . 

In this work, we study a particle condensation phe- 
nomenon in a driven pair exclusion process (PEP). The 
PEP was first introduced in Re f. (1211 a s a model for hub 
formation in evolving networks [13l - ll5l | . In the PEP, par- 
ticles hop under the so-called pair exclusion constraint, 
which will be explained later. With symmetric hopping, 
the system in the stationary state exhibits an intrigu- 
ing condensation state characterized by multiple meso- 
scopic condensates of size m con ~ TV 1 ' 2 and of number 
^Vcon ~ TV 1 / 2 , where N is the total number of particles 
with a logarithmic correction [12| ■ Those condensates 
are distributed randomly without any spatial correlation. 
We will investigate the effect of the hopping bias on the 
nature of mesoscopic condensation. 

The paper is organized as follows: In Sec. |TT] we in- 
troduce the PEP with and without a hopping bias. The 
condensation phase transition of the driven PEP has been 
investigated via numerical simulations, results of which 
are presented in Sec. IIII1 We summarize the paper in 
Sec. El 



II. DRIVEN PAIR EXCLUSION PROCESS 

There are M particles on a d-dimensional hypcrcubic 
lattice of N = L d sites under a periodic boundary con- 
dition. The occupation number at site i is denoted as 
rrii. Let us assume that there are M/2 distinct particle 
species and that each species has two elements. Pair ex- 
clusion means that no pair of particles of the same species 
is not allowed to stay at a same site. Such an interaction 
appears naturally in evolving networks |15| . 

The particle hopping dynamics is given as follows: At 
each time step, (i) we select a source site i at random 
from among N sites, (ii) When rrii > 0, we select one 
particles from among rm particles and attempt a particle 
hopping with the probability 



u(rrii) 



1 



/(I + 6). 



(2) 



(hi) There are 2d possible target sites {i ± ik\k = 
1, • • • ,d}, where e/t denotes the unit vector in the fcth 
direction. The target site is chosen from among forward 
sites {i + ik\k = 1, • • • ,d} with probability p or from 
among backward sites {i — e-k\k = 1, • • • , d} with proba- 
bility q = 1 — p. (iv) The hopping attempt is accepted 
only if it does not violate pair exclusion. The hopping 
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FIG. 1: The driven pair exclusion process in one (a) and two 
dimensions (b). Particle species are distinguished with filled 
patterns in (a). Pair exclusion forbids the particle represented 
by an empty circle to hop to the right. 



is symmetric when p = q = 1/2 while it is biased to the 
forward (backward) direction when p > q (p < q). The 
particle dynamics is illustrated in Fig. [TJ 

Due to pair exclusion, a hopping attempt is accepted 
or rejected depending on a particle's species distribution. 
In Ref. [12|, the accepting probability was shown to be 
approximately given by 



v(m) 
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where m < il'/ is the occupation number at a target 
site. Using the approximation, we can map the PEP to 
a driven diffusive system whose particle hopping proba- 
bility from site i to j is given by 
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p (q) to the for- 



with the hopping probability p^ 
ward (backward) direction. 

When the hopping is symmetric (p = q), the model is 
solvable, and its stationary state probability distribution 
is given by a product form as in Eq. ([1]) [12|. In fact, the 
model has the same stationary state as the ZRP with the 
hopping rate function given by 

b m 
uzRp(rn) = u(m)/v(m - 1) oc H h — . (5) 

771 M 

Note that the additional factor -S accounts for pair ex- 
clusion. Such a factor is irrelevant when m = 0(1), but 
it plays a crucial role in the condensation phase. It sup- 
presses a macroscopic condensate, and multiple meso- 
scopic condensates of size 77i con ~ TV 1 / 2 and of number 
^Vcon ~ N 1 / 2 appear [12|, [IlHllI • Those mesoscopic con- 
densates are spatially uncorrelated because the probabil- 
ity distribution has a factorized product form. 

In this work, we investigate the driven PEP with p = 1 
and q = 0. The hopping bias invalidates the mapping of 
the PEP to the ZRP; hence, the stationary state proba- 
bility distribution is not given by a product form. One 
obvious question is whether macroscopic or mesoscopic 
condensation occurs or not. Another interesting question 
is about a spatial correlation. Since the stationary state 
is not given by a factorized form, a spatial correlation 
in the particle distribution exists. Numerical simulation 
results on these issues are presented in the following sec- 
tion. 



III. NUMERICAL RESULTS 

Wc have performed extensive Monte Carlo simulations 
for the driven PEP in d = 1 and d = 2 dimensions. Par- 
ticles are allowed to hop only in the preferred forward 
direction (p = 1 and q = 0) (see Fig. [T]). We adopt the 
jumping probability in Eq. ([2]) with 6 = 4. This par- 
ticular value of 6 is chosen because it allows mesoscopic 
condensation for the symmetric PEP [12| ■ We start with 
M particles distributed randomly on N = L d sites, and 
data are measured in the stationary state over a time in- 
terval T > 10 8 . The system sizes are up to L = 16000 
in ID and L = 100 in 2D. The condensation transition is 
examined with the occupation number distribution 
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P ( m ) = Yd \^2 S ( m i> m ) 



which is the probability of a site having to particles. 



FIG. 3: Power-law scaling of the size of a typical condensate 
(6) "icon (□), the number of condensates N con (o), and the total 
number of particles belonging to the condensates m to tai (°) 
in ID (filled symbols) and 2D (empty symbols). Parameter 
values are b — 4 and p = 4. Solid line are guides for the eye. 




FIG. 2: The occupation number distribution P(m) is shown 
at various system sizes N = L d for ID in (a) and 2D in (b). 
Parameter values are p = 4 and b — 4. The inset shows 
P(m) for systems of size N = 16000 in ID and N = 100 2 
in 2D when the particle density p = 1/4, 1/2, and 4, which 
are below, equal to, and above the critical density p c = 1/2, 
respectively. The symbols represent the critical distribution 
of the corresponding ZRP given in Eq. ([7| with b = 4. 

When the particle density is low, the distribution func- 
tion decays exponentially (see the insets in Fig. [2]) in to. 
Without pair exclusion, the PEP reduces to the ZRP, 
which does not show condensation at low particle density. 
Since pair exclusion suppresses condensation, naturally, 
condensation does not occur. 

When the particle density is so low that there is no 
condensate in the system, the effect of pair exclusion rep- 
resented by Eq. ([3]) can be negligible in the infinite size 
limit. Namely, the PEP becomes equivalent to the ZRP 
in the normal phase without condensates. We expects 
that the equivalence persists up to the critical density 
p c = 1/(6 — 2) at which the ZRP undergoes a condensa- 
tion transition J4|. In the insets of Fig. ©, wc compare 



the occupation number distribution at p = p c with the 
critical occupation number distribution of the ZRP |4| 



P c {m) ex 
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with T(x) = (x — 1)!, which scales as P c (m) ~ m~ h for 
large m. The numerical data arc in good agreement with 
the critical distribution with b = 4. This comparison 
shows that the driven PEP undergoes a condensation 
transition at the same threshold p c = 1/(6 — 2). 

When the density is high, a broad peak in P(rn) ap- 
pears. The peak position moves to the right as L in- 
creases, as shown in Fig. [2] The peak represents con- 
densates. The peak is broad, so we quantify the size of a 
typical condensate, m con , with the highest peak position. 
The dependence of the condensate size on the number of 
lattice sites is plotted in Fig. [3] It follows a power-law 
scaling as 



m con (A0 ~ N? 



(8) 



with /3 ~ 0.56 in ID and /3 ~ 0.62 in 2D. The condensate 
is not macroscopic but mesoscopic with < j3 < 1. 

We note that the spectral weight of the condensate 
peak is much larger than 1/7V, which implies that there 
are multiple condensates. We quantify the number of 
condensates, iV con , from the total spectral weight of the 
peak. There is a local minimum in P(m) separating the 
occupation number into two regions. We estimated 7Vc 0n 
as the total spectral weight beyond the local minimum 
multiplied by N. Numerical data in Fig. [3] show that it 
also follows a power-law scaling as 



N r , 



N a 



(9) 



with a ~ 0.51 in ID and a ~ 0.34 in 2D. The total 
number of particles belonging to the condensates is pro- 
portional to N. 
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FIG. 4: (a) Occupation number distribution function C(r) 
around the largest condensate in ID with 6 = 4 and p = 4. 
The inset shows a snapshot of an occupation number distribu- 
tion, (b) Scaling plot of C{r)/N^ against r/N s with /3 = 0.56 
and 8 = 0.30. 



such condensate clusters. In fact, the inset of Fig. |4ja) 
shows that there are a few condensate clusters. With 
small values of (a — S) ~ 0.2, the expected number of 
clusters is small, and we could not verify the power-law 
scaling numerically. 
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We have shown that the driven PEP exhibits meso- 
scopic condensation characterized by the scaling m con ~ 
N* 3 and N con ~ N a . It contrasts with the CMA and the 
TP in which condensation does not occur in the presence 
of a hopping bias. Although the driven PEP exhibits 
a similar type of mesoscopic condensation as the sym- 
metric PEP, significant differences exist. The symmetric 
PEP can be approximated as a driven diffusive system 
with the hopping probability given in Eq. (j4]) , which can 
be mapped to the ZRP. Since the stationary state of the 
ZRP is given by a factorized product form, no spatial 
correlation exists. Consequently, the ZRP in all dimen- 
sions has the same property [12| . The variation of a and 
j3 with respect to d in the driven PEP suggests that a 
spatial correlation does matter. 

The spatial correlation is clearly seen from the spatial 
distribution of the occupation number. When conden- 
sation occurs, one can locate the site %m at which the 
occupation number is maximum. Then, one can measure 
the mean occupation number C(r) at site j = im + r dis- 
placed from %m by r. The distribution function is plotted 
for ID at b = 4 and p = 4 in Fig. g] (a). The occupation 
number is high near r = and then decays to a constant 
value at large r. This shows that mesoscopic condensates 
are bound to each other to form a cluster. We speculate 
that the clustering may originate from the jamming of 
driven condensates due to pair exclusion. 

The occupation number distribution turns out to fol- 
low a scaling form 



C{r) = N^T(r/N 5 ) 



(10) 



with /3 = 0.56 and 6 = 0.30 (see Fig. g^b)). The scaling 
function decays exponentially. There are N con ~ N a 
condensates, so one may expect that 5 — a. However, the 
numerical analysis in Fig. H[b) yields a value of 5 = 0.30, 
which is smaller than a = 0.51. The discrepancy may be 
explained by the assumption that there are N c i ~ N a ~ 5 



FIG. 5: Scaling of the occupation number along the directions 
(a) parallel and (b) perpendicular to the bias in 2D with b — 4 
and p = 4. 

We also study the clustering of mesoscopic condensates 
in 2D. Measuring the occupation number distribution 
C(r) at all lattice sites j = %m + t with r = xe\ + ye2, we 
found that all condensates form a single cluster. Inter- 
estingly, the condensate cluster has an elongated shape. 
Presented in Fig. [5] are numerical data C||(x), the occu- 
pation number at sites j = %m + x(e\ + e-i) along the 
direction parallel to the bias, and Cj_(cc), the occupation 
number at sites j = %m + x(ei — I2) along the direction 
perpendicular to the bias. They are fitted well to the 
scaling form 



C l{>± {x) = N^Q lx (x/N s ^) . 



(11) 



The exponent (3 = 0.62 is the same as the one for the 
condensate size. On the other hand, the two exponents 
describing the characteristic width of the cluster have dif- 
ferent values of Sn ~ 0.25 and 5± ~ 0.1. The total num- 
ber of mesoscopic condensates inside the cluster scales as 
N s \\+ S ± m ^0.35^ wm ch j s c l se to the previous estimate 
N a with a = 0.34. 

We add a few remarks on the shape of the occupation 
number distribution function. We observe that the dis- 
tribution is symmetric as C(r) = C(—r) in ID. In 2D, the 
distribution in the transverse direction is trivially sym- 
metric as C±(x) = Cj_(— x). However, it is asymmetric 
in the longitudinal direction as Cm (a:) 7^ Cn(—x). The 
origin of the asymmetry is not understood yet. 

From the spatial distribution of the occupation num- 
ber, we conclude that there is a clustering of condensates 
in the driven PEP. Condensate clustering was observed 
in interacting particle systems [T^, [2(| , where the clus- 
tering was caused by a particle attraction. Clustering 
in the driven PEP has a different origin. Clustering docs 



not occur in the unbiased PEP where mesoscopic conden- 
sates are distributed randomly. When there is a hopping 
bias, there is congestion due to the pair exclusion. The 
jamming leads to clustering. 



IV. SUMMARY 

We have investigated particle condensation in the 
driven PEP in ID and 2D. The PEP with symmetric 
hopping displays mesoscopic condensation, which is char- 
acterized by uncorrelated condensates of number N con ~ 
N a and of size m con ~ TV 3 with a = /3 = 1/2 in all di- 
mensions [12| ■ The driven PEP also exhibits mesoscopic 
condensation. However, a fundamental difference exists. 
The driven PEP is characterized by a spatial correla- 



tion that is manifested in the <i-dcpcndcnt exponents and 
the clustering of condensates. The clustering is a conse- 
quence of jamming caused by the combined effect of the 
hopping bias and the pair exclusion. The shape of the 
condensate cluster follows the scaling form of Eqs. (p~0|) 
and (fTTj) . It is interesting to note that the cluster is 
anisotropic in 2D. Our results show that the spatial cor- 
relation leads to rich behaviors of condensation phenom- 
ena, which need to be investigated thoroughly in the fu- 
ture. 
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